Multiloop integrals in dimensional regularization made simple 
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Scattering amplitudes at loop level can be expressed in terms of Feynman integrals. The latter 
satisfy partial differential equations in the kinematical variables. We argue that a good choice of 
basis for (multi-)loop integrals can lead to significant simplifications of the differential equations, and 
propose criteria for finding an optimal basis. This builds on experience obtained in supersymmetric 
field theories that can be applied successfully to generic quantum field theory integrals. It involves 
studying leading singularities and explicit integral representations. When the differential equations 
are cast into canonical form, their solution becomes elementary. The class of functions involved is 
easily identified, and the solution can be written down to any desired order in e within dimensional 
regularization. Results obtained in this way are particularly simple and compact. In this letter, we 
outline the general ideas of the method and apply them to a two-loop example. 



INTRODUCTION AND SUMMARY 

Scattering amplitudes are of fundamental interest in 
quantum field theory, as they connect theory to (collider) 
experiments. A theoretical challenge is to describe events 
with many jets, as well as to increase the precision of the 
predictions by going beyond leading and next-to-leading 
order accuracy, see [lj for a recent discussion. 

Computing scattering amplitudes usually consists of 
two steps. The first is to find an expression in terms of 
(Feynman) loop integrals, and the second is to evaluate 
the latter. In supersymmetric theories, the first step is 
essentially solved by virtue of (generalized) unitarity or 
recursion relations [2441 ; where the idea is to reconstruct 
the loop integrand from its analytic behavior, in particu- 
lar factorization on propagator poles. In developing these 
ideas, the maximally supersymmetric Yang-Mills theory 
was instrumental. What can we learn from recent ad- 
vances in this theory at loop level? 

Although many impressive perturbative results are 
available, unfortunately the methods used to obtain them 
cannot be applied directly to generic scattering ampli- 
tudes. The reason is that most of them relied in some 
aspect on the duality between scattering amplitudes and 
Wilson loops specific to that theory. Nonetheless, these 
advances suggest to us that there should also be eas- 
ier ways to directly carry out the Feynman integrations. 
Progress in this direction commenced with a better un- 
derstanding of loop integrands. It was observed that new 
representations for the latter |4|,|5| lead to simple analytic 
answers [g, H, as a consequence of differential equations 
they satisfy [8[. Moreover, mathematics for iterated inte- 
grals [9l4ll| furthered the understanding of the transcen- 
dental functions involved. 

In this letter, we apply these ideas to the evaluation of 
generic loop integrals in dimensional regularization, and 
combine them with state-of-the art techniques of QCD. 
Our method applies to planar or non-planar, massless or 
massive integrals equally. 

As we review shortly, the calculation of an arbitrary 
loop integral can always be related to the calculation of 



a finite set of master integrals. The question we wish 
to address here is how to choose a 'good' set of master 
integrals. 

One key idea is that we would like loop integrals to 
have simple properties under the action of differential op- 
erators. In order to define what is meant by simple, let us 
introduce the concept of degree of trans cendentality T(f) 
of a function /, which applies to a large class of iterated 
integrals [9J,llO(. T(f) is defined as the number of iterated 
integrals needed to define the function /, e.g. T(log) = 1, 
T(Li„) = n, etc. We also have T(fif 2 ) = T(/i) + T(/ 2 ). 
Constants obtained at special values are also assigned 
transcendentality, e.g. T(( n ) = n. Algebraic factors 
have degree zero. Functions that we will be interested 
in have uniform (degree of) transcendentality, i.e. if / is 
a sum of terms, all summands have the same degree. 

Moreover, we call such functions pure if their degree of 
transcendentality is lowered by taking a derivative, i.e. 
T{d f) = T(f) — 1. This implies that the transcendental 
functions in / cannot be multiplied by algebraic coeffi- 
cients, which would otherwise be 'seen' by the differential 
operators. 

There are several guiding principles that can help to 
find such integrals. Although to the best of our knowl- 
edge there is no general proof, it has been observed that 
integrals having constant leading singularities )5lll2| have 
these properties. The latter are defined by analytically 
continuing the momenta to complex values, and replac- 
ing the integration over space-time by contour integrals 
around poles of the integrand. Another way the prop- 
erties discussed above can be made manifest is when an 
appropriate 'd-log' representation is available [13(, where 
the integrand is written as a logarithmic differential form. 
This approach works particularly well for Wilson loop in- 
tegrals. Finally, we have also found explicit representa- 
tions based on Feynman parameters to be useful. 

One might think that few Feynman integrals have such 
nice properties. In fact, as we argue here, quite the op- 
posite is the case. This claim is supported by many ex- 
amples in the literature. As a specific case in point, the 
above ideas were employed in [14j to present massless pla- 



nar and non-planar form factor integrals in a basis where 
each integral has uniform transcendentality. 

Having discussed these general ideas, let us return to 
the issue of the reduction to master integrals and the 
calculation of the latter. 

Feynman integrals can be classified according to their 
topology, starting with the integrals where the maximal 
number of propagators is present. The propagators, la- 
belled by i, are raised to powers Oj. Subtopologies, where 
certain propagators are absent, are obtained by setting 
the corresponding indices a, to zero. 

For each topology, there is a set of integration-by-parts 
(IBP) identities [l5| that relates integrals with differ- 
ent values of the at. These equations follow from the 
Poincare invariance of the integrals, which is preserved 
in dimensional regularization. They are linear in the in- 
tegrals, with the coefficients being rational functions of 
the kinematical invariants and the space-time dimension. 
In this way, one can relate an integral with general inte- 
ger powers to a finite set of master, or basis integrals. In 
practice, a set of basis integrals can be found straightfor- 
wardly by using various public computer codes [lfj. For 



a recent review and more details and references, see [17 1. 

Having reduced the general problem in this way, one 
would then like to compute the basis integrals. From the 
discussion above it should not be surprising that we will 
use the method of differential equations [18- 21] . The idea 
is to differentiate w.r.t. the kinematical invariants. This 
can be implemented on the Feynman integrals by defining 
appropriate derivatives w.r.t. the momenta (respecting 
momentum conservation and on-shell conditions.) The 
r.h.s. of such an equation involves integrals of the same 
topology, but with different powers Oj. The latter in- 
tegrals can then be re-expressed via the IBP relations 
in terms of the chosen basis integrals. So in general we 
obtain a set of linear, first-order partial differential equa- 
tions for the basis integrals. 

Denoting the kinematical variables by Xi, the set of N 
basis integrals by fi, and working in D = 4 — 2e dimen- 
sions, this set of equations takes the form 



C^n/( e ) %n) — A m [e,x n )j(e,x n ) , 



(1) 



where d m — -*^—, and each A m is an N x N matrix. 



The matrices have to satisfy the integrability conditions 
(except possibly for special singular values of the x m ), 



<J n A m — u m A n + [A n , A m \ — , 



(2) 



where [A, B] := AB - BA. 

In order to completely specify the solution, one has 
to provide a boundary condition. This can in general be 
done by studying physical limits of the scattering process 
under consideration. 

In practice, one would like to solve eq. ([T]) in a Laurent 
expansion around e = 22 1. 



Under a change of basis / 
transform as 



Bf, the matrices A T 



A r , 



B~ 1 A m B 



B-\d m B) 



(3) 



Note that here B can be in principle any N x N matrix, 
where each entry is a function of e and of the kinematical 
variables x n . 

Equation ([T]) can simplify considerably when a good 
choice of basis is made. Our conjecture is that an op- 
timal choice of integral basis fi can be reached, where 
the integration of the system of differential equations be- 
comes trivial, in the sense explained below. 

One could imagine several simplified versions of eq. 
dU). We studied various cases of practical interest, and 
based on that evidence we propose that one can trans- 
form eq. ([IJ into the following form, 



d m f(e,x n ) =eA m (x n )f(e,x n ). 



(4) 



We remark that in this case, the integrality condition @ 
becomes 



<JnA m — O m A n — , [Ant A m \ — (J . 



(5) 



For the discussion of the properties of the solution it is 
convenient to combine eqs. (J4J and write them in differ- 
ential form, 

df(e,x n ) = edA(x n )f(e,x„) (6) 

We may also note that one can formally solve eq. (j6|) in 
terms of a path-ordered exponential, 

f=Pe*fc dA f(e = 0), (7) 

where the integration contour C connects the base point 
(representing the boundary condition) to x n . In other 
words, the perturbative solution in e is given by iterated 
integrals, where the entries of dA determine the integra- 
tion kernels. 

We stress that once a form (H|), or equivalently © is 
reached, then the problem of solving for the master inte- 
grals fi in the e expansion is essentially solved. 

The form (|6|) of the equations can also make the tran- 
scendentality properties of the solution manifest. In or- 
der to see this, let us introduce one new concept. In di- 
mensional regularization, it is customary to assign degree 
of transcendentality —1 to e. In this way, one can dis- 
cuss the transcendentality properties of functions appear- 
ing in the Laurent expansion of e-dependent expressions. 
For example, in x e = 1 + eloga; + 0(e 2 ) each summand 
has the same degree of transcendentality, namely zero. 
Then, if d A in ([6]) is a logarithmic one- form, it is clear 
that the answer will have uniform transcendentality, to 
all perturbative orders in e. 

Along the same lines, one can immediately determine 
from eq. ^ which class of functions the solution will be 



expressed in. This is best discussed using the notion of 
symbol of a transcendental function [9h11|. The entries 
of A in eq. ^ determine the alphabet of the symbols of 
the solution, again to all orders in e. 

In principle, starting with a random basis of master 
integrals /, one could attempt to find an appropriate set 
of functions B in eq. ([3j> in order to reach the canonical 
form (|4|). However, this seems like a formidable task in 
general, and for that reason the criteria outlined above 
are very useful in practice. 

In the remainder of this letter, we present a nontrivial 
example of this method. A more detailed discussion, as 
well as further applications, will be given elsewhere. 
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FIG. 1. Integral basis corresponding to /i,.../4 (first line) 
and /s,.../s (second line), up to overall factors. Fat dots 
indicate doubled propagators, and the dotted line an inverse 
propagator. The incoming momenta are labelled in a clock- 
wise order, starting with p\ in the lower left corner. 



Planar two-loop master integrals for 2 — > 2 scattering 



with the constant matrices 



We consider the planar double ladder integrals [21|,|23j. 
One can see via IBP that these constitute all loop inte- 
grals for virtual corrections to massless 2 — » 2 scattering, 
in any gauge theory. We introduce the notation 



h u 



:= e 2e7E 



A ^n(%))-, (8) 



{i-K D / 2 f 



■m—1 



with the propagator P(q) := l/(— q 2 ), and the set of pos- 
sible momenta q m , corresponding to TO = 1, ... 9, respec- 
tively, is ki,ki +pi,ki +pi2,ki +Pi23,k 2 ,k 2 +Pi2,k 2 + 
P123, ki - k 2 , where p 12 = p\ +p 2 and p 12Z = p x +p 2 +p 3 . 
We have p 2 = and J2 i — 1 Pi =0. The results depend on 
the Mandclstam variables s = 2pi ■ p 2 and t — 2p 2 ■ p 3 . 

There are 8 master integrals for this problem. We 
choose the following basis, see Fig. [TJ 



h = — e (~ s ) c t /o, 2, 0,0,0,0, 0,1,2 j 
h = € ( — s) e A), 0,2, 0,1,0,0, 0,2 j 
h = € {— s) e Io, 1,0,0, 1,0,1,0, 2 j 
/ 4 = -e a (-*) a + ae 7 2 , ,: 



.1,0.2,0.1,0.0 > 



\l+2e. 



h — € (—s) 1 " 1 "^ 6 4/l,l,l,0,0,0,0,l,2, 
fe = ~ e ( — s ) c \ s + t) ^0,1,1,0,1,0,0,1,1 

f _ A / \2+2e . t 

J7 — ~ £ \ — s ) t-'l, 1,1, 0,1, 0,1, 1,1 1 

.18 — ~ e \ s ) ^1,1,1,0,1,-1,1,1,1 • 



(9) 
(10) 

(11) 
(12) 

(13) 
(14) 
(15) 
(16) 



Here the normalization was chosen such that all fi admit 
a Taylor expansion in e. Moreover, they depend on s and 
t through the dimensionless variable x — t/s only. For 
this choice of basis, we hnd 



d x f = e 



/• 



(17) 
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and 



b = 
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Equation (|17j) is a simple instance of the Knizhnik- 
Zamolodchikov equations 24 1 25j. The three singular 
points {0, —1, 00} for x correspond to the physical limits 
s = 0, u = — s — t = and t = 0, respectively. 

Taking into account that the leading term in the e ex- 
pansion must be a constant, it follows from eq. (|17p that 
the result at any order in e can be written as a linear com- 
bination of harmonic polylogarithms [26| of argument x, 
with indices drawn from the set {0,-1}. In particular, 
the symbol alphabet in this case is {x, 1 + x}. Finally, the 
requirement that the planar integrals be finite at x = —1 
and real-valued for x positive turns out to fix all except 
two boundary constants. The latter can be related to 
the trivial propagator- type integrals f 2 and /4, which are 
known in closed form. This completely solves this family 
of Feynman integrals, to all orders in e. We see that all 
basis elements fi have uniform degree of transcendental- 
ity, to all orders in e. One may verify agreement with ref. 
[23j, to the order in e computed there. 



DISCUSSION AND OUTLOOK 

In this letter, we argued that a good basis choice for 
master integrals can significantly simplify the differential 
cqu; lions they satisfy. We motivated and discussed guid- 



ing principles tor choosing good master integrals, based [4 

on their transcendentality properties. The latter can be 
discussed through their leading singularities. We remark 
that the basis choice is not unique. 

We provided a non-trivial example at the two-loop level 
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anylorder in e. All master integrals have uniform degree 
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leads to (TJ|. We would like to stress that beyond the 
example given here, we found this method to apply to 
many further cases of practical interest. In particular, 
we expect applications to previously unknown integrals 
involving top quarks or Higgs particles, or to Bhabha 
scattering, to name a few examples. Preliminary results 
also show that the method can be applied successfully to 
integrals in heavy quark effective theory. 



hi more complicated multi-leg processes, or processes 



H9 
[10 

[11 



[12 
[13 



[14 



involving masses, the appropriate set of integral functions rjg 



may not yet be known. We anticipate that our method 



will be a convenient way of solving this problem, and lead 1 16 
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